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The electronic structure of nanocylinder without and with a small perturbation is investigated 
with the help of calculation of the local density of states. A continuum gauge field-theory model 
is used for this purpose. In this model, Dirac equation is solved on a curved surface. The local 
density of states is calculated from its solution. The case of 2 heptagonal defects is considered. This 
paper is an extension of our previous work [J] where one heptagonal and one pentagonal defects in 
hexagonal graphene network were compared. The metaUization for the perturbed cylinder structure 
is found. 
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I. INTRODUCTION 



The carbon nanostructures play a key role in constructing nanoscale devices like quantum wires, nonlinear electronic 
elements, transistors, molecular memory devices or electron field emitters. Their molecules are variously-shaped 
geometrical forms its surface is composed of disclinated hexagonal carbon lattice. The main structure of this kind is 
graphene - the carbon lattice plane from which all other kinds are derived. 

The most famous form is fuUerene - a material composed of molecules which have the form of a soccer ball 
0. Other kinds are nanocones, nanotubes, nanotoroids, nano cylinders, nanoribbons etc. The various forms of 
the nanostructures are ensured by the topological defects in the graphene which are most often presented by 
the pentagons and heptagons in the hexagonal plane lattice. In the closest vicinity of these defects, the positive 
resp. the negative curvature arises for the case of the pentagons or the heptagons, respectively. Generally, 
for the n-sided polygon, n < 6 corresponds to the positive curvature and n > 6 corresponds to the negative 
curvature. This fact coincides with the choice of the form of the defects in the particular cases. There is no 
heptagonal defect in the fuUerene, but a lot of these kinds of defects can be found in many open forms of nanostruc- 
tures. But we can find them also in some of the closed forms like nanotoroids or more complicated, folded forms 
of nanotubes. Most often, the heptagons appear in pairs with pentagons in the connecting parts of the folded forms 

Because of the applications, the research of the electronic properties of the carbon nanostructures is important. 
One of the main characteristics is the local density of states (LDoS). In the presented model coming from 
the effective-mass theory, knowledge of the solution of the corresponding Dirac equation is necessary for the 
calculation This solution is represented by the wave- function and to find it, we have to know the geometry 
of the molecular surface. From the mentioned facts follows that the chosen geometry can be only an approxi- 
mation of the complicated real situation: for example, the spherical geometry of fuUerene is not suitable for the 
I description of the closest vicinity of the defects but it correctly describes the properties of the whole molecule. 
As discovered in fE] for the case of nanocones, the most suitable geometry for the description of the close vicinity 
of the defects is the hyperboloidal geometry. Very often, for a given geometry, the number of possible defects is limited. 

The solutions for spherical, conical and 2-fold-hyperboloidal cases were found in 0-0] ■ In [H, we used the presented 
model for calculation of the electronic properties of the structures with the geometry of the 1-fold hyperboloid. 
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The aim was to describe the electronic properties in the vicinity of the locally negative curvature of an arbitrary 
nanoparticle. This restriction did not enable us to do the calculations for the case of more than 1 defect. In this 
paper, we present a model describing the electronic properties of a simple nanocylinder and a curved nanocylinder 
including 2 heptagons at the opposite sides of the surface. The hyperboloidal geometry is used again. Because a 
nanocylinder is an opened nanotube, comparison with the case of the capped nanotube could be performed (see e.g. 
[3] for this purpose). 

This paper is organized as follows; the second section describes the computational formalism. The third section 
summarizes the basic properties of nanotubes and derived nanostructures. The fourth section researches electronic 
properties of cylinder without and with a defect. In the fifth section, a small review about the graphene nanoribbons 
and their properties is given and calculation of its LDoS is performed. In Conclusion, the obtained results for cases 
of cylinder and inifinitely long nanoribbon are compared and discussed and a brief review about the production of 
the nanocylinders is introduced. For the cylinder, the normalization constants are computed in the Appendix A and 
zero modes for defect-free and perturbed case are computed in the Appendices B and C, respectively. 



II. COMPUTATIONAL FORMALISM 

To research the electronic properties, we have to solve the Dirac equation in (2+1) dimensions. It has the form 

ia"e^ [d^ + ia^ - ia^]iP = E^, (1) 
where means the partial derivation according to the parameter, i.e. = gfTr- 

In this equation, besides the energy E, the particular constituents have the following sense: 
ct", a = 1, 2, denote the Pauli matrices. 

The zweibein e^,/i = z^ip stands for incorporating fermions on the curved 2D surface and it has to yield the 
same values of observed quantities for different choices related by the local SO (2) rotations: 



ea^^e'^^A^ep, Af e 50(2). (2) 

For this purpose, a covariantly-constant local gauge field is incorporated Q: 

d^.^: - r^,e^ + iu:,)y, = 0, (3) 

where 

VP ^laP- (^ + ^^^] (4) 

2^ \dxP dx^ dx^ ) ^ ' 

is tli6 Levi-CivitcL connection coming from tlie metrics g^j^y (see below). Then is called the spin connection. 
The constituent 

n,^^u;f[ac,,a0] (5) 

denotes the spin connection in the spinor representation. Its components are 

fiz = 0, fl<p = iujcjs, (6) 

where 
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The sense of the metric coefRcients will be explained below. 

The wave function ■(/;, the so-called bispinor, is composed of two parts: 



B 



each corresponding to different sublattices of the curved graphene sheet. The gauge field arises from spin rotation 
invariance for atoms of different sublattices A and B in the Brillouin zone jlOl and the gauge field is connected 
with the chiral vector (n, m) [III [13 : 
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= N/A, = -^{2m + n). (9) 



The metric 17^,^ of the 2D surface comes from following parametrisation, with the help of two parameters z, 
ip: 

{z, tf)-^!^ = {x{z, (f), y{z, 'p),z), (10) 

where 

0<Z<OO, 0<<y9<27r. (11) 

The 4 components of the metric are defined as: 

9^.. = d^ltd.lt. (12) 
For the rotationally symmetric case which will be researched, the non-diagonal components of the metric are 

ff^v = g^z = 0- (13) 



If we write the wave function in the form 



and substituting (ITil) into ([1]) we obtain 



j = 0,±l,... (14) 



-_u — £jV, ■ V — hu, (15) 



where 



j=j + 1/2 ~a^~ c^. (16) 



Each of the solutions u, v consists of two linearly independent components such that 

u{E, z) = Ci{E)ui{E, z) + C2{E)u2{E, z), 



(17) 



\ ( '^zUl JUi \ C2 I dzU2 JU2 \ 

E y^/gZ yfW^ ) E y^/gTz ^fW^ j 
where for a concrete value of E^ the functions Ci{E), C2{E) stand for satisfying the normalization condition 



(18) 



27r / {\u(E,z)\'^ + \v{E,z)\'^)Az = l. 



(19) 



FIG. 1: Nanocylinder with a small perturbation 



For a given zq, the LDoS is defined as 

LDoS{E) = \u{E, zo)|' + \v{E, zo)|^ (20) 



III. CARBON NANOTUBES AND NANOSTRUCTURED CYLINDERS 

Nanostructured cylinders consist of rolled-up graphene sheet from which the nanotubes arise by adding cups to the 
open edges. The properties of carbon nanotubes are presented in They can be characterized by the chiral vector 
defined as 

= not + m(i^ — {n, m), (21) 

where n,m d and at, a| are the unit vectors of the chosen coordinate system of the original graphene plane [l3j . 
The direction and length of coincide with the circumference of the nanotube. 

According to the symmetrical properties, we distinct 3 kinds of the carbon nanotubes: 

• armchair (ac); Zt — [n,n), 

• zig-zag (zz): Zt = (n,0), 

• chiral: otherwise. 

One can easily see that ac and zz are the only variants with the mirror symmetry. 

Furthermore, the carbon nanotubes can be divided into metals or semimetals depending on whether or not n — m 
is a multiple of 3. Evidently, all kinds of ac nanotubes are metals. 

In this paper, we are concerned with the case of cylinder without and with a small perturbation, caused by 2 
heptagonal defects (see Fig. [1]). 
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IV. LOCAL DENSITY OF STATES 



The following calculations of LDoS suppose long cylinders with large diameters. The z axis is identical with the 
axis of the cylinder and the values of z coordinate go from —z^ax to z^ax- It means that for z = 0, the results must 
be similar to the solution for the simple graphene and for z = ±Zmax, the results should correspond to the edge states 
as stated in But as we will see, our method of calculation is not accurate enough to distinct the edge states 

for the ac and the zz variants, respectively. The reason is that the Dirac equation involves information about the 
geometrical structure but not about the space coordinates of the individual atoms. 

The calculations will be first done for the defect-free cylinder and then for the perturbed cylinder. 



A. Defect-free cylinder 

The defect-free cylinder can be described with the help of the parametrization 

l^{z,(p) = {acoa(p,asm(p, z), (22) 
where a is the radius. The metric coefficients will be then 

gzz = 1, = (23) 

and 

detg^^=g = a^ (24) 
Because the metric coefiicients are constant, i.e. — 0, it follows for the coefficients of the Levi-Civita connection 

' ^ ox^ ' 

and the parameter w related to the spin connection 

r^. = 0, 0. = i. (25) 



From this follows 



so 



^ -u = Ev, -^-^, = Eu, (26) 



where 



The solution has the form 



'9zz \j9w \j9zz \j9w 



d,u - -M = Ev, -d^v --v = Eu, (27) 

a a 



J = j + 0.5 + i(2TO + n). (28) 



^(z) = Ci exp(Q!z) -I- C2 exp(— az), (29) 



v{z) = ^ ( ~ ~ I exp(Q;z) — -^ \ a+ — \ exp(— az), (30) 



E \ a I ' ' ' E \ a 

where 



a = a{E) = W 4 - ^2 (31) 



a2 
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FIG. 2: LDoS as a function of _B G ( — 1, 1) and z = 98 for defect-free cylinder 




FIG. 3: LDoS as a function of -E G (—1, 1) and z G (0, 100) for defect-free cylinder 

The procedure of the calculation of the normalization constants Ci (E) , C2 (E) is described in the Appendix A. An 
interesting task is the case of the zero modes which corresponds to the states of electrons at the Fermi level. They 
are derived in the Appendix B. 

In Fig. [2]the LDoS is shown as a function of E in the interval (—1, 1). In Fig. ^LDoS is plotted as a function of 2 
variables, E and z. It is case of the armchair configuration. As aforementioned, long cylinders with large diameters are 
supposed, so we choose high values, i. e. m = 20, n = 20, a = 14, z^ax = 100 and j = 0. Because of the exponential 
character of LDoS, its values are negligible for a long range of z values, so we choose a high z in the 2D plot of Fig. 
[21 i. e. 2; = 98. For the same reason, similar values will be chosen for the perturbed case. 
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B. Perturbed cylinder 



In the case of a small perturbation, we have 

'^{z,ip) = ^oVT+"A0^cos(/7, a-\/T+"A^sin</5, 2;^ , (32) 

where A is a positive real parameter, A << 1. For A = 0, we get the defect- free cylinder discussed in previous 
chapter. 
Then 

= 1 + ^^^J ~ 1 + a^AV, g^^ = a^l + Az^) (33) 

and 

g = g^^ = a^{l + Az^). (34) 
The only nonzero coefficients of the Levi-Civita connection will be 

Tl.^a^A^z, r^^~-a2A2, T^^ = T^, ~ a^A^ (35) 

and the oj parameter 



w~ ^(l-aA^^;). (36) 



Then we solve the system of equations 



, -'-^^=Ev, - , -l^^=Eu, (37) 

VTT^W^ a^l + Az^ v^T^W^ a^i + Az^ 

where 

~ 1 N 

j=j + l/2+-{2m + n)--, (38) 

i.e. we consider the perturbation to be caused by N heptagonal defects. For small A and neglecting the second order 
of A, it can be simplified as 

d^u - ^ (^1 - " = -^^^ ^^^^) ^ " ^^^^ 

The solution is 

u{z) = CaiD,, (e(2)) + Ca2^.. m^)), (40) 



v{z) = %i (d.D^Mz)) - I^^(l - Ia^.^)) + %a (a.I?.,(^^(z)) - 



i^-(!iM(i_lAV) , (41) 

CI A I 



where 



vi = i- 



V2 = —i- 



. o^A - Aa^E^ - 4zaVSj + 4j2 



8a^/Aj 



(42) 



,(.) = (-A)VM^/J+^|.l, 



(43) 
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FIG. 4: LDoS as a function oi E e (-1, 1) and 2 = 98 for perturbed cylinder with A = 0.05 (left) and A = 0.1 (right) 





FIG. 5: LDoS as a function of S £ (-1, 1) and z £ (0, 100) for perturbed cylinder with A = 0.05 (left) and A = 0.1 (right) 

Dv{C) being the parabolic cylinder function [l^. The functions Cai = Cai{E), Ca2 = C/\2{E) will be calculated in 
the same way as for the defect-free cylinder, i. e. from the normalization condition 



{\u{E,z)f + \v{E,z)f)dzdip^ATT / {\uiE,z)f + \v{E,z)f)dz ^1 



(44) 



In Fig. m a comparison of the 2D plots of LDoS is made for different values of A and similarly in Fig. [5l where 
the 3-D plots are compared. The number of defects N = 2 and the other values are same as in the defect-free case. 
On the whole, we can conclude from the comparison that LDoS is slightly increasing for higher A. Although this is 
only the first order approximation, we will see that the real solution of p7p behaves in a similar manner. The case 
of zero modes appropriate to the states at the Fermi level is calculated in the Appendix C. 

Comparing the calculated expressions (|B1[) . (jC2l) for zero modes corresponding to the defect-free and perturbed 
case, respectively, we see that it has a character of an exponential function related to a polynomial in z its order 
depends on the chosen precision. This means that the simplest case of the first order polynomial corresponds to the 
defect-free case, the higher orders of the polynomial correspond to the perturbed case. For higher z, the function 
uo{z) corresponding to the sublattice A is increasing, whereas the function vq{z) corresponding to the sublattice B 
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is decreasing. 



V. CASE OF NANORIBBONS 



For large diameters, the presented results for the defect-free cylinder should correspond with the density of states 
for nanoribbons (ac or zz). The reason is that because of a very small curvature, a nanocylinder with a large diameter 
we can take as an approximation of an infinitely long nanoribbon (see e.g. [l6j). 

Let us give a small review about the structure and properties of the nanoribbons. Graphene nanoribbons are 
materials, theirs molecules have the form of thin, rectangular pieces cut from the graphene layer, so their structure 
is analogical. Similarly as in the case of graphene nanotubes, ac and zz forms are distinguished, theirs difference 
consists in the form of the edges. Denoting TV — 1 the number of hexagonal unit cells by which the nanoribbon width 
is constituted, metallic and semiconducting forms can be marked again: a concrete form of nanoribbon is metallic, 
if iV + 1 is a multiple of 3, so this condition differs from the condition mentioned in the third chapter for carbon 
nanotubes. 

To calculate the energy spectrum, the tight-binding approximation is used. Deriving the form of Hamiltonian for 
the given variant (ac or zz), we get the set of equations of motion. For ac, it has the form p7| 

where m = 1, . . . , denotes the position of the corresponding A or B atom in the nanoribbon, ipm.A or ^m,B, resp. 
is its wave function and k — ^n, n — 0, ±1, . . . , ±(-^ — 1), {Ly is the number of the hexagonal unit cells in the 
longitudinal direction), is the longitudinal wave number. With the help of the boundary conditions and respecting 
the fact that the determinant of the system (j45|) should be equal to zero, we get that the energy spectrum contains 
the values 




£: = ±Jl + 2e^cos - +£2, (46) 



where Er — "2 cos jf^T^, r = 1, 2, . . . iV. Using this procedure for zz, the analogous energy spectrum appears, but 
the formulas are a little bit modificatcd. In [18], similar results are obtained for the carbon nanotubes. 

Now the LDoS for ac nanoribbons will be calculated and the results will be compared with the above stated 
calculations for the defect-free nanocylinder. The corresponding formulae is 

LDoS^E.,,^'^^^. ,47) 

where 



D{E) = lim 2Im / dk — , (48) 

5{E — E{k)) is the Dirac delta function and k is the wavenumber. Let us stress that the integration over k is used 
here instead of the integration over z. Both kinds of integration are related by the Fourier transform, as follows from 
For different k, we see the plots of LDoS in Fig. [51 Evidently, excluding the case k = 0, the results more or less 
correspond to the solution given by the formulas (P^ . (150)) . 

In the possible solutions for zz are divided into the cases of the extended and the localized states. The form 
of the corresponding wavefunction for the localized states is 

Since for large z we can write 



sinh z ^ exp z, (50) 
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FIG. 6: LDoS as a function of S G (^1) 1) for fc = (left), k — ^ (middle) and k = n (right) for ac nanoribbons 



this form corresponds to the solution for the defect- free case (|29l) . (1501) . So, we see an important coincidence with the 



results found for nanoribbons in [17| . The results presented in this paper are related to the case of ac configuration. 



but as stressed in the beginning of the fourth section, our method of calculation does not distinguish between ac and 
zz and the plots for case of zz would be similar. 



Because of the discrete spectrum of carbon nanoribbons and from the fact that it can be sometimes understood as 
an approximation of the nanocylinders a question arises about a possibility of calculation the low energy electronic 
states for the case of small defects. Similar task was done in jl9] for the case of the spheroidal geometry. This 
could be possible in case of using different methods for calculation of the LDoS. In this paper, the resulting energy 
spectrum can be understood as continuous, so there is not any possibility of shifting the energy levels due to the 
defects. 



VI. CONCLUSION 



Comparing the plots of LDoS in Fig. [6] with the plots in Fig. [2l we see that both Figures are similar but not 
completely same. The reason is in difference between the used procedures: for the purpose of Fig. [6l the energy 
spectrum is first calculated in (j46|) and then the corresponding LDoS is found according to ((47|) . but the procedure 
used in chapter IV is opposite, i.e. first the LDoS is calculated and then the energy spectrum is estimated from the 
plots (as in [1]). Because such a procedure is only a rough approximation, it does not enable us to find the localized 
states for zero energy as in (H71) . So, only the method described in the previous chapter enables us to distinguish 



metals and semimetals. Furthermore, according to [17[, the condition for nanoribbons to be metal is different from 



the condition for nanotubes, so the results could never be exactly same. 

The presented results for perturbed cylinder are only the first order approximation of the real solutions which are 
presented in Figs. [7]and|Sl Because of the complicated form of ([571) . the solutions were acquired in a numerical way. 
The analytical solutions we can get in a more precise form, if we increase the order of the parameter A in p9p . It 
follows from the plots on logarithmic scale in Figs. [71 [51 (and also from Fig. 21 for the first order solution) that for 
E = 0, a localized state with a small amplitude appears. From this we see the evidence that the metallic properties 
of the nanocylinder are increasing by adding a perturbation into the structure. 

To conclude, the LDoS was researched for nanocylinders without and with a defect. We found that for the 
defect-free case, the electronic properties can sometimes correspond to the case of infinitely long nanoribbons. For 
the case of perturbation, the metallic properties are more manifested. 

The generation of the nanocylinders is in a close correspondence with the generation of the nanotubes. The most 
exploited method is the chemical vapor deposition (CVD) which is based on the production of the required material 
from the surface of a substrate by volatile precursors. Usually it is realized as a thermal reaction called catalytic 
growth process, when the growth of carbon nanotubes is catalyzed in a heated flow furnace by a carbon surface 
diffusion [lO, [2lj . This process is activated by iron and nickel nanoparticles with the help of acetylene and hydrogen 
as the precursors. The nanocylinders are produced in this process together with the nanotubes or they arise from the 
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FIG. 7: LDoS as a function of iJ G (—1, 1) and z £ (0, 100) for perturbed cylinder with A = 0.05 on linear (left) and logarithmic 
(right) scale (numerical solution) 




FIG. 8: LDoS as a function of E £ ( — 1, 1) and z G (0, 100) for perturbed cylinder with A — 0.1 on linear (left) and logarithmic 
(right) scale (numerical solution) 

nanotubes with the help of other thermal or chemical processes. The surface of the nanotubes is usually um-eactive, 
but the change of its structure can be provoked by catalytically active atoms of rhenium inserted into the nanotubes 
[2^ . For the case of the perturbed nanocylinder, as presented in this paper, its structure can be determined by the 
presence of two heptagonal defects on the surface, but this structure sometimes appears without defects and it is only 
given by different lengths of the bonds between the atoms. The calculations performed in this paper are valid anyway. 

Appendix A: Normalization constants of defect-free cylinder 



We will suppose that for the solution holds 

u{z) = u{—z), v{z) = v{—z). 



(Al) 
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Then the normalization condition gives 



{\uY + \v\')Azdip^2 j {\u\' + \v\')2T:dz 





Sttj 



a a j 



{Cl{E){^- - a(S))[exp(2«(^)z„,,) - 1] - Cl{E){l + a(S))[exp(-2a(i;)z„„,) - 1]}+ 



aa{E)E' 



167rC\iE)C2{E)z^a.^l- 



From (I29l).(l30l) and (IM 



(A2) 



LDoS{E, z) ^ \u{E, z)|2 + \v{E, : 



- + 



+Ci{E) - - a{E) — 2exp(2a 

\ a / aE^ \ 2azmax 



+ CiEf I + a{E) -j^ (2 eM-Mm + ^^P(2«w) 1 
\ a / aE^ \ 2az„,n:r 



(A3) 



Now, if we suppose Ci{E) e TZ for aU E e (-1, 1), it foUows from (|A2 



aa{E)E^ 



\ 47rj(i - a(^))(exp(2a(S)z„,,) - 1) 



(A4) 



So, for the given choice of the parameters, we find the dependence / = f{E) and we can choose 
Ci{E) = min{f{E),E E (—1, 1)}. Usually, we choose Ci{E) ~ C2{E), but this choice is not necessary. 



Appendix B: Zero modes for defect-free case 



The solution of the zero modes has the form 



uo{z) = Cu exp I j I voiz) = Cy exp | -^z 



In accordance with the previous procedure, 

4n J (\uo{zr + \vo{z)\')dz^^lcl 



exp I ^rnax 

\ a 



ct 



exp I ZfYidx 
\ a 



and 



LDoS{z) = C'l 



2j_ 

exp I ^maa; 

\ a 



j/(27ra) 



exp -fz, 



(Bl) 



1 (B2) 



(B3) 
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Then 



Cu< 



j7(27ra) 



\ exp (^^ 

f2 TfldX j 



(B4) 



Appendix C: Zero modes for perturbed case 



The zero modes for the case of perturbation will be calculated from the equations (|39| without the r.h.s, i.e. 



d,uo - Ml - ^Az2 j uo = 0, -d,vo - ^1 - ^Az^ ) wo - 0, 



(CI) 



with the solution 



uo{z) = Cau exp 



vo{z) = Cav exp 



(C2) 



We require 



An J (|Mo(z)r + |«oWr)dz = 47r J <( exp 





Cav exp 



dz 



= 47r / < C^„exp ( ^z j exp(-fc^) +C'i„exp I -^z j exp((523) i dz = 1, 



(C3) 



where we denoted ^ — 6- In case Sz^^^ << 1, an approximate expansion can be done 



exp((5z^) ^l + Sz^ 



(C4) 



and the last integral equals 



47r 



Ci„exp I j (1 - Sz^) + Ci„ exp | -^z J (1 + 5z^) \ dz 



aTT ^2 J / 2jz„ 



3(5(a - jz„a:E)(a^ - ajZmax+fz'^^^)+f{4: - (5z^„^) 



3a3(5 - 4/ 



^^^M^^H a 



3S{a + ]z„,ax){a' + ajz„,ax + j'^^aJ + r (4 + fc™ax) - 3a^(S - 4j^ ^ = 1. (C5) 



If (Jzj^j^^ << 1 is not true, we do the expansion 



exp(i5z^) = 1 + (5z^ 



3\n 



2! 



+ ... + 



(C6) 



where ^'^^ } << 1, substitute it into the integral and make an approximate calculation. Calculations of LDoS, Cau 
and Cav are then done in the usual way. 
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